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L'invention concerne le domaine des appareils destines aux therapies par 
hyperthermic locale. L'invention concerne egalement un precede de mise en ceuvre de 
tels appareils. 

Les therapies par hyperthermic locale consistent a chauffer, localement, une 
zone cible d'un tissu biologique. Lorsque ce type de therapie est utilise dans le cadre de 
la therapie genique, la chaleur peut etre, par exemple, utilisee pour son action sur un 
promoteur thermosensible. On peut egalement utiliser la chaleur pour necroser des 
tissus biologiques et proceder a l'ablation de tumeurs. 

Aussi, les therapies par hyperthermic locale offrent-elles de nombreux 
avantages. Ces avantages sont aussi bien qualitatifs qu'economiques. Du point de vue 
qualitatif, elles offrent, par exemple, un fort potentiel pour le controle de traitements tels 
que les therapies geniques, le depot local de medicaments, l'ablation de tumeurs, etc. 
Du point de vue economique, elles sont compatibles avec un traitement ambulatoire des 
malades, elles permettent de reduire les durees d'hospitalisation, etc. 

Dans les therapies par hyperthermic, la chaleur peut etre, par exemple, apportee 
par un laser, des micro-ondes ou des ondes radiofrequences, des ultrasons focalises, etc. 
D'une maniere generate, les therapies par hyperthermie locale permettent des 
interventions medicales dont la nature invasive est reduite au minimum. Mais parmi les 
types d'energie precites, les ultrasons focalises sont particulierement interessants 
puisquMls-permettent-de-chauffer-la-zoncde-Jbcalisation, de maniere non invasive, 
profondement dans un corps biologique, sans chauffer de maniere significative les tissus 
voisins de la zone de focalisation. 

Dans tous les cas, la temperature de la zone cible et de son environnement 
immediat, pendant le traitement, doit etre controlee precisement et de maniere continue, 
alors que l'apport en energie est localise et rapide (de l'ordre de quelques secondes). A 
cette fin, on peut installer des sondes de temperature dans la zone cible et son 
environnement immediat. Mais on peut egalement utiliser l'Imagerie par Resonance 
Magnetique (ERM). L'IRM permet, en effet, d'obtenir une cartographie precise des 
distributions de temperatures ainsi que des informations anatomiques detaillees. L'IRM 
permet en outre un controle non invasif de la temperature. 

On connait deja des dispositifs pour le controle de la temperature, pendant les 



traitements par des ultrasons focalises, bases sur la thermometrie par resonance 
magnetique. De tels dispositifs sont en particulier decrits dans les documents suivants: 
"Control system for an MRI compatible intracavitary ultrasons array for thermal 
treatment of prostate disease", Smith NB et aL, Proceedings of the annual meeting of the 
5 International Society of Magnetic Resonance in Medicine, 1999, p. 672 et "Real time 
control of focused ultrasound heating based on rapid MR thermometry", Vimeux FC et 
aL, Invest. Radiol. 1999, 34, P- 190-193. Dans les dispositifs decrits dans ces 
documents, le retrocontrole de la chaleur apportee par les ultrasons focalises, grace aux 
cartographies obtenues par IRM, est de type PDD (acronyme de 1' expression anglo- 

10 saxonne Proportional Integral and Derivative). De plus, avec ces dispositifs, le controle 
de la chaleur fournie au tissu, repose sur la prise en compte d'une temperature mesuree 
dans la zone de focalisation des ultrasons uniquement, ou correspondant a une moyenne 
obtenue a partir de la distribution spatiale de la temperature, dans la zone cartographiee. 
La figure 1 represente revolution temporelle de la temperature moyennee de la 

15 zone de focalisation, traitee grace au dispositif decrit dans le premier de ces documents. 
Sur cette figure, la temperature evolue jusqu'a un plateau correspondant a la 
temperature que Ton souhaite atteindre dans la zone de focalisation. On peut constater 
que la temperature souhaitee dans la zone de focalisation n'est atteinte qu'apres une 
duree de l'ordre de 30 minutes. 

20 La figure 2 represente revolution temporelle de la temperature moyennee dans la 

zone de focalisation , traitee grace au dispositif decrit dans le deuxieme des documents 

cites ci-dessus. On peut constater que la temperature souhaitee dans la zone de 
focalisation est atteinte en moins de 2 minutes. Par contre, on observe des variations de 
la temperature souhaitee, de plus ou moins 4° C. 

25 Un but de Tinvention est de proposer un ensemble pour le traitement thermique 

d'une zone cible d f un tissu biologique permettant a la fois d'obtenir rapidement la 
temperature souhaitee dans la zone cible et de maintenir et controler la temperature dans 
cette zone cible avec une precision accrue, par rapport a ce qui etait possible avec les 
techniques de Part anterieur. 

30 Ce but est atteint, selon 1'invention, grace a un ensemble pour le traitement 

thermique d'une zone cible d'un tissu biologique, comprenant 

- des moyens generateurs d'energie pour fournir localement de Tenergie dans la zone 



cible ; 

- des moyens pour mesurer et enregistrer la temperature dans la zone cible; 

- une unite de controle comprenant des moyens pour determiner, a partir de la 
temperature mesuree dans la zone cible, la quantite d'energie devant etre foumie a la 
zone cible, et des moyens pour commander aux moyens generateurs de l'energie, de 
delivrer cette valeur de la puissance ; 

caracterise par le fait que l'unite de controle comprend en outre des moyens de 
traitement numerique, point par point, de la distribution spatiale de la temperature dans 
la zone cible et son environnement. 

L'ensemble de traitement thermique conforme a la presente invention, prend en 
compte la distribution spatiale reelle de la temperature dans la zone cible, mais aussi 
dans l'environnement de cette zone. C'est a dire qu'il prend en compte et traite, point 
par point, cette distribution spatiale. Contrairement aux ensembles de traitement 
thermique de Part anterieur, la distribut ion spatiale de la temperatureest utilisee pour en 
deduire des gradients de temperature et non pas de simples moyennes. Ceci permet 
d'estimer, avec plus de precision, quelle doit etre la quantite d'energie a deposer et done 
de parvenir plus rapidement a la temperature souhaitee et de maintenir la temperature du 
tissu biologique avec une plus grande stabilite. 

Avantageusement, l'unite de controle de l'ensemble de traitement thermique 
selon l'invention, comprend en outre des moyens pour estimer les pertes locales en 
— energie-thermique-a-partir-dlune-estimation_d e la conduction thermique et de la 



distribution spatiale de la temperature dans la zone cible et son environnement. En effet, 
^information fournie par la valeur des gradients de temperature, ainsi que la prise en 
compte d'une estimation des pertes locales de chaleur, permettent non seulement de 
connaitre la facon dont le tissu biologique traite a reagi a la chaleur deja appliquee a 
celui-ci, mais permettent en outre, grace a la prevision sur la facon dont le tissu 
biologique va reagir a la chaleur. Ceci permet egalement de faire evoluer la temperature 
du tissu traite thermiquement, plus rapidement vers la temperature souhaitee et de 
maintenir la temperature du tissu biologique avec une plus grande stabilite. 

Avantageusement, les moyens generateurs d'energie de l'ensemble de traitement 
thermique selon l'invention, emettent des ultrasons focalises. En effet, les ultrasons 
focalises permettent de fournir de la chaleur, dans une zone tres localisee, de maniere 



non invasive, meme si cette zone est situee profondement dans un corps humain ou 
animal. En outre la focalisation permet de ne pas chauffer de fafon significative, les 
tissus avoisinant la zone de tissu biologique traitee. 

Avantageusement, les moyens pour mesurer et enregistrer la distribution spatiale 
5 de la temperature de P ensemble de traitement thermique selon Pinvention, comprennent 
un appareil d'imagerie par resonance magnetique. En effet, PIRM permet une mesure 
de la temperature, non invasive, precise, et bien resolue en de nombreux points de la 
zone cartographies Les donnees recueillies par IRM sont, en outre, aisement traitees 
numeriquement. 

10 Avantageusement, Tensemble de taitement thermique selon Pinvention comprend 

des moyens devaluation de la distribution spatiale, dans la zone cible et son 

environnement, de l'energie fournie a la zone cible. 

Selon un autre aspect, Pinvention est un procede de regulation d'un ensemble de 

traitement thermique d'une zone cible d'un tissu biologique, comprenant 1'etape 
15 consistant a deposer localement de l'energie dans la zone cible, 

caracterise par le fait qu'il comprend en outre les etapes consistant a 

- evaluer les gradients de temperature dans la zone cible et son environnement ; et 

- en deduire l'energie a deposer dans la zone cible pour atteindre la temperature 
souhaitee. 

20 Avantageusement, ce procede comprend alors en outre Pet ape consistant a 

estimer les pertes locales en energie, dans la zone cible et son environnement. 

Avantageusement, ce procede comprend en outre Petape consistant a evaluer la 
distribution spatiale, dans la zone cible et son environnement, de l'energie fournie a la 
zone cible. 

25 D'autres aspects, buts et avantages de l'invention apparaitront plus clairement a 

la lecture de la description detaillee qui suit. 

L'invention sera aussi mieux comprise a l'aide des dessins joints sur lesquels: 

- la figure 1 represente revolution temporelle de la temperature d'une zone cible, 
lorsque celle-ci est traitee en utilisant un ensemble de traitement thermique de Tart 

30 anterieur ; 

- la figure 2 represente 1'evolution temporelle de la temperature d'une zone cible traitee 
par un autre ensemble de traitement thermique de Tart anterieur ; 
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- la figure 3 represente schematiquement l'ensemble de traitement thermique conforme 
a l'invention ; 

- la figure 4 est un organigramme du precede de regulation conforme a la presente 
invention ; 

5 - la figure 5 represente 1'evolution du Laplacien en fonction du temps, au cours du 
traitement thermique par l'ensemble conforme a la presente invention, correspondant 
a une experience in vitro decrite ci-dessous ; 

- la figure 6 presente les resultats d'une serie d'experiences visant a etudier l'influence 
des erreurs d'estimation sur les pertes d'energie, representees par la diffusivite de la 

10 chaleur, et le coefficient d'absorption de 1' energie ultrasonore ; 

- la figure 7 represente 1'evolution temporelle de la temperature maximum, mesuree au 
cours du traitement thermique d'un tissu biologique, par l'ensemble conforme a la 
presente invention, correspondant a la meme experience in vitro que celle de la 
figure 5 ; 

- la figure 8 represente 1'evolution temporelle de l'amplitude du signal emis par le 
generateur de l'ensemble conforme a la presente invention, au cours du traitement 
thermique d'un tissu biologique, correspondant a meme experience in vitro que celle 
des figures 5 et 7 ; 

- la figure 9 represente 1'evolution temporelle de la temperature maximum, au cours du 
traitement thermique, d'un tissu biologique, par l'ensemble conforme a la presente 
Inventi on, cones pondant-a-une-autre-experieneew>»^//w^avec-trois-palierS-de_ 
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temperature, et decrite ci-dessous ; et 

- la figure 10 represente 1'evolution temporelle de la temperature maximum, mesuree 
au cours du traitement thermique par l'ensemble conforme a la presente invention, 

25 correspondant a une experience in vivo decrite ci-dessous. 

Un des modes de realisation de l'invention, est decrit ci-dessous de maniere detaillee. 
A titre d'exemple, ce mode de realisation de l'invention correspond a un ensemble de 
traitement par hyperthermic locale par ultrasons focalises, controles par IRM. 
Comme represente sur la figure 3, un tel ensemble comprend : 
30 - des moyens generateurs d'energie 100 ; 

- des moyens de cartographie 200 ; 



- une unite de controle 300 ; et 

- un porte-echantillon 400 pour le tissu biologique 410 a traiter. 

Dans le mode de realisation de Finvention decrit ici, les moyens generateurs 
d'energie 100 sont composes d'un transducteur 110, d'un generateur de signal 
5 sinusoidal 120, d'un amplificateur 130 et d'un convertisseur 140, reliant le generateur 
de signal sinusoidal 120, a l'unite de controle 300. 

Le transducteur 110 opere a 1,45 MHz. Un tel type de transducteur 1 10 est, par 
exemple, commercialise par la societe Speciality Engineering Associates® (Soquel, 
Californie). Son diametre et sa distance focale sont de 38 mm et 25 mm respectivement. 
10 Le generateur de signal sinusoidal 120 est, par exemple, du type FG110, 

commercialise par la societe Yokogawa® (Tokyo, Japon). 

L'amplificateur 130 est, par exemple, du type KMP 170F, commercialise par la 
societe Kalmus® (Bothell, Washington). Cet amplifiacteur 130 a un gain en puissance 
de 58 dB. 

15 Le convertisseur 140 est, par exemple un convertisseur BEEE488 serie, 

commercialise par la societe I. O. Tech® (Cleveland, Ohio). 

Les moyens de cartographie 200, permettent de mesurer et enregistrer la 
distribution spatiale de la temperature. lis comprennent, par exemple, uri appareil IRM 
de type Bruker Biospec commercialise par la societe Bruker® (Ettlingen, Allemagne). 

20 Cet appareil utilise un aimant de 4,7 T qui est equipe d'un insert de 120 mm de 

diametre, qui genere des gradients du champ magnetique (la valeur maximum du 

gradient est de 0,193 T/m). 

L'unite de controle 300 comprend en particulier une station de travail 310 de 
type Alpha PW 500a MHz, commercialism par la societe Digital®. 

25 L'unite de controle 300 comporte egalement des moyens devaluation et de 

traitement numerique de la distribution spatiale de la temperature 320, des moyens de 
determination de la valeur de la puissance 330 devant etre fournie a la zone cible, des 
moyens d'estimation des pertes locales de Tenergie thermique 340 et des moyens de 
commande 350 des moyens generateurs d'energie 100. Les moyens de commande 350 

30 indiquent aux moyens generateurs d'energie 100 de delivrer la valeur de la puissance 
fournie par les moyens de determination de la valeur de la puissance 330. 

Le porte-echantillon 400 comprend un support 420 de rat, en plexiglas®. Ce 
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support 420 contient le transducteur 1 10 et une bobine de surface (non representee sur la 
figure 3). Un tel porte echantillon 400 a deja ete decrit dans les documents « Fast lipid 
supressed MR temperature mapping with echo-shifted gradient echo imaging and 
spectral-spatial excitation », de Zwart JA et al., 1999, Magn ResMed, 42, p.53-59 ; et 

5 « On the feasibility of MRI-guided focused ultrasound for local induction of gene 
expression », Madio DP et al., 1998, J. Magn. Res. Imaging.I, 8, p. 101-104. Le support 
420 est place dans un tube de plexiglas® qui est partiellement rempli d'eau. Le 
transducteur 110 est positionne de maniere a ce que le point de focalisation 460 des 
ultrasons soit situe approximativement a 10mm en profondeur, dans le tissu biologique 

10 410. Pendant les mesures in vitro, une sonde 430 de temperature est inseree dans le tissu 
biologique 410 constitue d'un morceau de viande fraiche, de maniere a avoir une 
reference pour la temperature. Cette sonde 430 est par exemple un thermocouple du 
type Digi-Sence DualLog, commercialise par la societe Cole-Parmer Instrument Co.® 
(Vernon Hill, Illinois). 

15 La preparation des echantillons pour les experiences in vitro et in vivo est 

realisees de la maniere suivante. On prend des rats males de la race Wistar de 325 a 500 
g. Ceux-ci sont anesthesies en associant 1% en volume d'halothane avec un melange 
constitue de 7 volumes d'oxyde nitreux (N 2 0) pour 3 volumes d'oxygene, selon un 
protocole approuve. Afin d'ameliorer la penetration du faisceau d'ultrasons dans le tissu 

20 biologique 410, la cuisse du rat (qui se trouve entre transducteur et point focal) est 

epilee, *" ntilisant un produit prevu a cet effet et connu de l'homme du metier. Pendant 

les mesures in vivo, la temperature endorectal des rats est enregistree. La temperature 
corporelle des rats est maintenue a 35°C, par immersion du corps des rats dans un bain 
dont la temperature est regulee a cette fin. Apres le traitement thermique les rats sont 

25 sacrifies. 

Les moyens de cartographie 200 sont utilises en mettant en oeuvre une sequence 
Echo gradient, avec les parametres suivants : TR = 50 ms, TE = 15 ms, la taille de la 
matrice = 64x63, trois lignes de l'espace des k par TR, FOV = 64 mm x 63 mm, ou TR 
est le temps de repetition, TE est le temps d'echo et FOV est le champ de vue (FOV est 
30 Pacronyme de l'expression anglosaxone "Field of View"). Les donnees sont obtenues 
par IRM, a partir d'une coupe de 2 mm d'epaisseur, perpendiculaire au transducteur a 
ultrasons et comprenant le point focal des ultrasons focalises. La resolution temporelle 
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des cartographies obtenues par ERM est de 1,05 s. La resolution spatiale des 
cartographies obtenues par IRM est de 1x1x2 mm 3 . Les cartographies de la temperature 
mesuree par resonance magnetique sont obtenues a partir des mesures du deplacement 
de la frequence du proton de l'eau. Le choix de la resonance du proton de l'eau, pour 
5 faire ces mesures, repose sur le fait que la relation entre la frequence de resonance du 
proton de l'eau et la temperature est, en premiere approximation, independante de la 
composition du tissu biologique 410. Le deplacement de la frequence du proton de l'eau 
en fonction de la temperature est en outre lineaire et cette linearite n'est pas affectee par 
les modifications du tissu biologique 410 induites par la chaleur (Ishihara Y et al., 

10 Magn. Res. Med., 1995, 34, p. 814-823 ; Peters RD et al., Magn. Res. Med., 1998, 40, p. 
454-459). La dependance en temperature de la frequence de resonance du proton est de 
0,0094 ppm-K" 1 (« Fast magnetic-resonance temperature imaging » ;de Zwart JA et al, 
1996, J. Mag. Res. B, j_12, p. 86-90 ; « Fast lipid supressed MR temperature mapping 
with echo-shifted gradient echo imaging and spectral-spatial excitation », de Zwart JA 

15 et al., 1999, Magn. Res. Med, 42, p.53-59). 

Les signaux de resonance magnetique provenant des lipides constituent une 
source significative d'erreurs dans les cartographies de temperature calculees, puisque 
les frequences de resonance des protons des lipides ne dependent pas de la temperature. 
On supprime done les signaux de resonance magnetique provenant des lipides, en 

20 utilisant une excitation selective de l'eau, de la maniere decrite dans le document « Fast 
lipid supressed MR temperature mapping with echo-shifted gradient echo imaging and 
spectral-spatial excitation », de Zwart JA et al., 1999, Magn. Res.Med, 42, p.53-59. 

L'utilisation des moyens devaluation et de traitement numerique de la 
distribution spatiale de la temperature 320, a egalement deja ete decrit dans le document 

25 «Fast lipid supressed MR temperature mapping with echo-shifted gradient echo 
imaging and spectral-spatial excitation », de Zwart JA et al., 1999, Magn. Res. Med, 42, 
p.53-59. 

Un exemple de mise en oeuvre du procede de regulation de l'ensemble de 
traitement thermique, conforme a la presente invention est decrit de maniere detaillee 
30 ci-dessous. 

Sur la figure 4, est represents schematiquement, un organigramme de ce mode 
particulier de mise en oeuvre du procede selon T invention. 



Conformement a cet exemple de mise en oeuvre, le procede selon 1' invention 
comporte : 



- une etape 1 d'estimation du coefficient de diffusion de la chaleur (on) et du 
coefficient d'absortion des ultrasons focalises (a 2 ) dans le tissus biologique 410 ; 

- une etape 2 de definition par I'utilisateur, d'un profil de revolution temporelle 
de la temperature souhaitee ; 

- une etape 3 d'acquisition d'une image IRM ; 

- une etape 4 de calcul de la distribution spatiale de la phase au point focal 460 et 

dans son environnement ; 

- une etape 5 d'etablissement de la distribution spatiale de la temperature au 

point focal 460 et son environnement ; 

- une etape 6 devaluation des gradients de temperature du point focal 460 et 
dans son environnement ; 

- une etape 7 de determination de la nouvelle puissance devant etre delivree par 

le generateur 120 ; et 

- une etape 8 de changement du taux d'energie deposee par le generateur 120. 
Les etapes 3 a 8 sont reprises en boucle, afm d'atteindre et de parcourir le profil 

de revolution temporelle de la temperature souhaitee, defini a l'etape 2. 

Pour un element donne du transducteur 110, la puissance electrique P(t) 
tTmsmise-a-un-echanti^ 



de la puissance 320. Sa valeur peut done etre modifiee en direct par l'unite de controle 
300: Elle est obtenue sur la base de l'equation : 

m = -7=r4 ^ - «,0O- vr_(o +*•[©«)- • m) 

Cette equation est etablie en se fondant sur les considerations suivantes. 

Le point focal est defini par r= (0,0,0). Involution temporelle et la 
temperature maximum au point focal correspond alors a r msic (/)= 7(0, 0,0,/). Notons 
par 0(r) le profil predetermine de revolution temporelle souhaitee de la temperature 
maximum T^if). Comme indique ci-dessus, ce profil est defini avant le debut de 
chaque experience. Par exemple, ce profil ©(/) comporte une montee de 10° C, pendant 
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100 s. Cette etape de montee suit revolution d'une moitie de la periode de la fonction 
cosinus. Elle est suivie par une periode ou la temperature est constante (10° C au-dessus 
de la valeur de depart), pendant 250 s. 

La temperature maximum 3^(0= ^(P 9 0 9 0 9 i) peut etre controlee seulement au 
5 point focal. En effet, la geometrie du transducteur 110 et la distribution spatiale de 
l'indice de refraction dans le tissu biologique 410 determine le champ acoustique. En 
consequence, revolution de la temperature, ailleurs que dans la zone de focalisation, fait 
intervenir des fonctions dependantes de la coordonnee d'espace r et de la temperature 
T. Le champ de la puissance acoustique p(r ), le tenseur de la diffiisivite de la chaleur 
10 a,(r,7) dans le tissu biologique 410 et le coefficient d'absorption des ultrasons 
focalises a 2 (r,T), sont alors lies par la relation : 

^|^= ^ryVTir^y^T) p(ryP(t). [la] 

Lorsque la diffiisivite est isotrope et varie lentement spacialement, Pequation la 
15 se simplifie pour donner l'equation lb, ou a,(r,r)= ^7>[a 1 (r,r)] est un champ scalaire 

d 2 

et V 2 est l'operateur Laplacien defini par —7 + —7 +~TT : 

at cy oz 

= a, (r, T). V*T(r,<) + a 2 (r, T\ p{r)- P(t) [lb] 

ot 

II faut noter que les fonctions a,(r,7) et a 2 (r 9 T) ne sont pas precisement 
connues au debut du chauffage. 
20 Le champ de la puissance acoustique d'un element spherique de transducteur a 

ultrasons focalises, correspond approximativement a une distribution gaussienne autour 
du point focal 460, avec un rayon d'attenuation a 6dB note Ro* Le temps de diffusion 

D2 

specifique r est defini par r = . Pour une longueur d'onde ultrasonore de 1 mm, 

Tordre de grandeur de rest de Tordre de 10 s. 
25 On se fixe les objectifs suivants : 

1) la temperature du point focal 460 doit atteindre la temperature souhaitee, aussi vite 
que possible (c'est-a-dire dans un temps du meme ordre de grandeur que r), ceci sans 
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oscillation, ni depassement de la valeur souhaitee ; 

2) une fois cette valeur souhaitee atteinte, la temperature doit rester constante pendant 
une periode predefinie par l"utilisateur. 

L'integrale du profil 0(0, l'integrale de revolution temporelle de la temperature 
maximum observee experimentalement T mac (tlet la difference entre ces deux integrates, 
sont respectivement : 

Q(t) = \e(f)df Pa] 

0 

o 

a(0 = Q(0- v(0 = j MO- T^ye [2c] 

0 

En utilisant ces expressions, liquation lb, donnant Involution de la temperature 
au point focal, peut etre exprimee en fonction de Aft): 

tm . Ml . afj^y v'T^ W - «,&.> />(«) P] 

at at 
our = (0,0,0) est omis et /<0,0,0) = 1 . 

Dans l'equation 3, le parametre que Ton souhaite controler directement est la 
puissance des ultrasons focalises P(t). Notons qu'une equation differentielle d'ordre 2 et 
lineaire en A(r), peut etre avantageusement utilisee par l'unite de controle 300, de 
-TMniereTmitlaire-a^n-systeme-d 



A(r) d'une telle equation tend asymptotiquement vers zero, et qu'il en est de meme pour 
sa derivee premiere. Si la derivee premiere de A(f) est egale a zero, 7^(0 presente un 
recouvrement avec le profil predetermine de 1'evolution temporelle de la temperature 
©(/). Ceci constitue l'idee fondamentale du precede de controle mis en ceuvre par l'unite 
de controle 300. Ainsi, nous pouvons reecrire l'equation 3 sous forme d'une equation 
differentielle lineaire de second ordre en Aj/X-de type : 



<PA d& a 1 r 4 i 

— -+a—~ + A = 0 l H J 

dt 2 dt 4 

De maniere a obtenir l'expression voulue de l'equation 4 a partir de l'equation 3, 
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P(t) est reecrit avec l f < 



expression suivante : 



P(/) = 



1 [ d®(t) 
cz 2 (T nax )[ dt 



-«l(^max) 




[5] 



L'equation 5 correspond a l'equation centrale utilisee pour calculer en direct le 
niveau de puissance des ultrasons focalises. A partir de la solution de l'equation 

5 differentielle du second ordre, correspondant a l'equation 4, on peut voir que le 
parametre a est relie au temps de reponse caracteristique t r de la boucle de regulation, 
par l'expression a=2l t r . On suppose que dans les equations 4 et 5, toutes les fonctions 
utilisees pour calculer la puissance P(t) sont precisement connues. On peut egalement 
verifier, comme le montre T equation 6 ci-dessous, que la temperature observee 

10 experimentalement T(i), tend asymptotiquement vers le profil 0(/) : 



Comme deja note plus haut, dans une experience reelle, les coefficients 
d'absorption des ultrasons a2 et le parametre de diffusion de la chaleur ainsi que Ieur 
dependance en temperature, sont inconnus. Ces parametres a, et a 2 dependent de la 
15 composition du tissu biologique 410, de processus physiologiques, telle que la 
perfusion, et de changements irreversibles intervenant pendant la procedure de 
chauffage, par exemple dans les procedures d'ablation. Ainsi, un systeme de regulation 
doit il etre tolerant vis a vis des erreurs sur les parametres a, et a 2 . 



Tequation 5, deux difficultes apparaissent : 

1) ?^ax(0 e t ^T^if) telles qu'elles sont obtenues a partir de la cartographie de la 
temperature issue de TIRM, sont affectees par le bruit. 

2) les valeurs de a, et cc 2 et leur dependance en temperature ne sont pas precisement 
25 connues, de meme leur sensibilite a la necrose par chauffage (par exemple dans les 

procedures d'ablation) et les parametres physiologiques telle que la perfusion, ne sont 
pas connus avec precision. 



®(0 - r{t) = ^(/) = [©Co)- no)} Ci-f-0- exp r f ■ 0 



[6] 



20 



Seul le profil ©(/) et sa derivee sont precisement connues. 

Aussi, lorsque la puissance ultrasonore est calculee en direct a partir de 



Toute erreur pouvant affecter et a 2 peut etre traitee comme une erreur de 
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parametre dans une boucle de controle, selon un modele lineaire. Ainsi, des estimations 
des valeurs -initiates dea> et a, seront choisies, puis utilisees pendant la procedure de 
chauffage, pour le calcul de la puissance ultrasonore, conformement a l'equation 5. 

Une analyse theorique de l'effet de l'erreur sur les parametres «, et a 2 dans 
l'equation 5 met en evidence les effets suivants : 

1) Une estimation erronee du parametre a, diminue l'efficacite de la boucle de 
regulation; il peut en effet s'ensuivre une possibility de depassement de la 
temperature souhaitee ou une valeur de la temperature determinee trop faible ; ceci a 
alors pour consequence une augmentation du temps de convergence. Quoiqu'il en 
soit, meme dans ces conditions, la temperature experimentale tend toujours 
asymptotiquement vers le profil predetermine de Involution temporelle de la 
temperature. 

2) Une estimation erronee du parametre a, conduit a un decalage constant, dans la 
zone ou le profil 0(0 est plat, des valeurs de la temperature, entre les valeurs de la 
temperature mesurees experimentalement et le profil 0(/). Ce decalage est 
proportionnel a la derivee premiere par rapport au temps, du Laplacien, multipliee 
par 1'erreur absolue sur a, , fois a 2 . Afin d'estimer cet effet, la derivee du Laplacien 
a ete determinee en utilisant une regression lineaire de la courbe representee sur la 
figure 5, entre 150 et 250 s. Sa valeur est de l'ordre de 0.01 K mrn^.s' 1 , ce qui 

20 conduit a un decalage en temperature approximativement egal a 0.1 °C. Ainsi, 

Veneur sur la temperature reelle ne devrait pas etfe~directemenm3eWeendu-fairde- 
la limitation de la precision des mesures thermometriques par resonance 
magnetique, a cause du bruit. 

L'influence de l'erreur sur les valeurs a x et a 2 sur l'efficacite du controle de la 
25 temperature par 1'ensemble de traitement thermique conforme a l'invention est etudiee 
experimentalement, sur un morceau de viande fhriche, ci-dessous. 

Selon I'etape 1, les parametres «, et Oj doivent etre estimes. Ceci est realise, 
dans un premier temps a partir d'une experience preliminaire avec une puissance 
d'ultrasons focalises constante. Le parametrea! est calcule a partir de la derivee par 
rapport au temps, de la temperature au point focal, divisee par la moyenne (sur cinq 
images IRM) du Laplacien. II est calcule immediatement apres l'extinction de 
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remission des ultrasons focalises et est exprime en mm 2 /s. Le parametre a 2 (la vitesse 
de depot de l'energie au point focal prenant en compte la distribution spatiale de la 
puissance des ultrasons focalises) est calcule a partir de la derivee, au temps initial, de la 
temperature du point focal, par rapport a la puissance des ultrasons focalises (lorsque les 
5 ultrasons focalises sont emis et que la diffusion est negligeable, voir P equation 1). II est 
exprime en K.s^.(mV)"2 La precision estimee pour a, et a 2 est meilleure que 10%, 
comme on peut le deduire d' experiences repetees. Les valeurs numeriques obtenues de 
cette fa9on peuvent etre directement utilisees dans P equation 5 pour calculer, a chaque 
fois qu'une nouvelle cartographie de la temperature est disponible, la valeur reelle de la 

10 puissance devant etre fournie par le generateur 120. 

Selon Tetape 2, un profil de revolution temporelle souhaitee de la temperature 
au point focal est defini avant le debut de chaque experience. Ce profil comporte une 
partie initiale croissante, correspondant a une moitie de la periode de la fonction 
cosinus, suivi d'une partie constante en temperature. La derivee premiere de la courbe 

15 correspondant a ce profil, est continue et peut etre calculee numeriquement par 
ordinateur. 

L'etape 3 est mise en oeuvre grace aux moyens de cartographie 200. 

Les etapes 4 a 6 sont mises en oeuvre grace aux moyens devaluation de 
traitement numerique de la distribution spatiale de la temperature 320. Les phases sont 
20 calculees a partir des signaux IRM, obtenus grace aux moyens de cartographie 200. Les 
changements de frequence de resonance deTeau sont calculees a partir de ces phases. 
Les changements de temperature sont calcules a partir de ces changements de 
frequence. 

La temperature maximum T max (t) et Tintegrale A(t) sont directement deduites des 
25 cartographie obtenus par IRM, respectivement au cours des etapes 5 et 6. L'integrale 
A(t) (Equation 2c) est calculee numeriquement a 1'aide de la station de travail 310. 

Au cours de Tetape 6, les pertes locales de l'energie thermique sont reevaluees a 
Taide des moyens d'estimation des pertes locales de Tenergie thermique 340. 

Selon Tetape 7 du procede conforme a la presente invention, Poperateur 
30 Laplacien V 2 T(r y t) est applique au traitement des cartographies de la temperature 
obtenues par IRM. La valeur de ce Laplacien au point focal est calcule, en utilisant la 
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methode des elements fins, en combinaison avec un filtre temporel de reduction du 
bruit. Ce filtre utilise un rapport de ponderation binomial egal a 1 :4 :6 :4 :1, sur cinq 
images. Cette etape 7 est mise en oeuvre avec les moyens de determination de la valeur 
de la puissance 330. 

Finalement, la puissance est calculee grace a 1'equation 5 et la valeur calculee est 
envoyee au convertisseur 140 (etape 8). L'etape 8 est mise en oeuvre avec les moyens de 
commande 350 des moyens generateurs d'energie 100. 

La duree totale de calcul pour traiter chaque cartographie de distribution spatiale 
de la temperature, c'est a dire chaque cycle de l'ensemble des etapes 3 a 8 decrites ci- 
dessus, est inferieure a 250ms. 

Afin d'analyser la tolerance vis-a-vis des erreurs dans les estimations initiales de 
a, et a 2 , plusieurs procedures de chauffage ont ete mises en oeuvre avec les parametres 
a, etaj variant sur une large gamme. Cette gamme va de 0 a 300 % d'une valeur 
predetermines pour a, , et de 40 a 250 %, pour a 2 . Cette valeur predeterminee est celle 
obtenue a partir des mesures preliminaires exposees ci-dessus. Des temps d'attente de 30 
minutes ont ete introduits entre les experiences successives dans le but d'atteindre les 
lignes de base en temperature, spatialement uniformes dans l'echantillon et identiques 
pour chaque experience. Neuf resultats representatifs sont reportes sur la figure 6. Ces 
resultats montrent que le systeme a une large tolerance vis-a-vis des erreurs sur les 
estimations des valeurs initiales de ^ et a r Quoi qu'il en soit, on peut noter que la 
l^ciende-controte 

instability est exacerbee lorsque est fortement surestime et a 2 est sous-estimee. Cet 
effet peut etre attribue au bruit experimental des mesures de la temperature par IRM. 
Lorsqu'un calcul affecte par le bruit, conduit a une surestimation de la valeur du 
Laplacien (les derivees secondes etant sensibles au bruit), la puissance des ultrasons 
focalises appliques augmente dans une proportion egale a a h sla 2y ou e est la 
surestimation du Laplacien. Une augmentation de la puissance des ultrasons focalises 
conduit a une forte augmentation du Laplacien dans le tissu biologique 410 et de ce fait, 
a une nouvelle augmentation de la puissance des ultrasons focalises. Cette reaction 
positive cessera apres un temps approximativement egal a 2/a, grace a une reaction 
negative dans la boucle de controle. Ceci explique la periodicite des instabilites dans ce 
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cas extreme. 

La force de la reaction negative sur la boucle de regulation correspond au 
parametre a. En effet, comme nous Tavons vu plus haut le parametre a, est egal a deux 
fois Tinverse du temps de reponse caracteristique t r (a=2/ t r ). La valeur de ce parametre 
5 est indiquee pour chaque exemple de la figure 6. En general, les valeurs de a de 0,1 s" 1 a 
0,2 s" 1 sont suffisantes pour atteindre des temps de montee en temperature similaires a 
ceux du profil ©(/), meme lorsque des valeurs fortement erronees de a, et a 2 sont 
utilisees. C'est seulement dans le cas extreme oil aj est fortement sous-estime et a 2 est 
surestime (voir figured, en bas a droite), que a doit etre augmentee jusqu'a la valeur de 

10 0,40 s" 1 pour obtenir un recouvrement avec le profil predetermine de revolution 
temporelle de la temperature. La valeur de a optimale, experimentalement trouvee (dans 
tous les cas, sauf les cas d'erreurs extremes sur a, et or 2 ), est de 0,2 s* 1 , le temps de 
reponse de la boucle de regulation correspondante etant de 10 s. Lorsque la force de la 
reaction negative est augmentee, une correction plus rapide des erreurs sur les 

15 parametres initiaux est obtenue, mais Tamplitude de la puissance des ultrasons focalises 
et des fluctuations en temperature autour de la valeur predeterminee est aussi 
augmentee. 

La mise en oeuvre et les performances de Pensemble de traitement selon 
Pinvention decrit ci-dessus, sont illustrees ci-dessous a I'aide de deux exemples. 
20 Exemple 1 : Utilisation de l'ensemble de traitement thermique, conforme a la 
presentejiwention,dans4ecadre^e4nesur^^>i^//w 



Selon cet exemple, un protocole de montee en temperature de 10°C est mis en 
oeuvre sur un echantillon de viande fraiche. La temperature initiale est egale a 15°C. 
Avec ce protocole, aucune modification irreversible du tissu biologique 410, resultant 

25 de ce profil devolution temporelle de la temperature, n'est attendue. La figure 7 
represente revolution de la temperature maximum en fonction du temps. Dans la partie 
plate de cette courbe, Televation moyenne de la temperature est de 9.97°C, avec un 
ecart type egal a 0.1 9°C. Cet ecart type doit etre compare a celui, egal a 0.1 8°C, qui est 
obtenu par les mesures de temperature efFectuees sans chauffage par ultrasons focalises 

30 (c'est a dire ce qui correspond a la ligne de base du bruit dans les mesures de 
temperature). Sur la figure 5, est reporte le Laplacien calcule en direct. L'attenuation 
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observee sur la partie constante en temperature correspond a la diminution des gradients 
en temperature autour du point focal. L'amplitude de la puissance appliquee en direct 
est reportee, quant a elle, sur la figure 8. A cause du bruit de mesure, la valeur calculee 
du Laplacien et l'amplitude de la puissance fournie par le generateur, presentent une 
fluctuation de l'ordre de 10%, approximativement. Ceci n'a qu'un faible effet sur la 
temperature resultante, car la frequence de fluctuation (c'est a dire l'inverse de la 
resolution temporelle de la cartographie par resonance magnetique) est beaucoup plus 
grande que l'inverse du temps de reponse (r) specifique au chauffage du tissu biologique 
410. 

La figure 9 montre la stabilite en temperature obtenue avec un profil a trois 
paliers (15, 25, 30 °C). L'ecart type est egal a 0.35°C, 0.36°C et 0.40°C, respectivement 
pour des montees en temperature a 15°C, 25°C et 30°C. Les resultats represents sur la 
figure 9 confirment la grande stabilite en temperature, sur une large gamme de montee 
en temperature, du systeme de regulation de l'ensemble de traitement thermique 
conforme a la presente invention. 

Exemple 2 : Utilisation de l'ensemble de traitement thermique, conforme a la 
presente invention, dans le cadre de mesures in vivo. 

En adoptant une procedure analogue a celle mise en oeuvre dans le cas de 
T exemple 1, des experiences ont ete menees in-vivo, sur une cuisse de rat. Les resultats 
correspondants sont montres sur la figure 10. La resolution temporelle est de 0.5s. La 
-t^mperature-moyennerentre-90-et~l-20s-apres-!e-deb 
valeur du profil a atteindre est egale a 55°C) avec un ecart type de 0.33°C. Les figures 
7, 9 et 10 montrent que Ton peut controler la temperature avec une precision proche de 
celle donnee par les mesures de temperature effectuees in vitro ou in vivo. 

II a ete decrit ci-dessus un mode de realisation de Tinvention correspondant a un 
ensemble de traitement par hyperthermic locale par ultrasons focalises, controles par 
IRM, mais l'invention couvre une gamme plus large d'ensembles de traitement 
thermique. Aussi, on comprendra que l'invention peut etre generalisee aux cas ou la 
chaleur est, par exemple, apportee par un laser, des micro-ondes ou des ondes 
radiofrequence, des ultrasons focalises, etc. On comprendra egalement que d'autres 
moyens de mesure de la temperature, peuvent etre utilises dans 1' ensemble de traitement 
thermique conforme a la presente invention, a la place de PIRM. 
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De meme revaluation et la traitement numerique de la distribution spatiale de la 
temperature ont ete decrits ci-dessus comme etant effectues a l'aide du Laplacien. 
D'autres moyens d'effectuer cette evaluation peuvent etre utilises sans sortir du cadre de 
T invention ; 
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REVENDICATIONS 

1. Ensemble pour le traitement thermique d'une zone cible d'un tissu biologique 
(410), comprenant 

- des moyens generateurs d'energie (100) pour fournir localement de l'energie dans la 
zone cible ; 

- des moyens (200)pour mesurer et enregistrer la temperature dans la zone cible; 

- une unite de controle (300) comprenant des moyens (330) pour determiner, a partir de 
la temperature mesuree dans la zone cible, la quantite d'energie devant etre fournie a la 
zone cible, et des moyens pour commander (350) aux moyens generateurs de l'energie 
(100), de delivrer cette valeur de la puissance ; 

caracterise par le fait que 1'unite de controle (300) comprend en outre des moyens de 
traitement numerique (320), point par point, de la distribution spatiale de la temperature 
dans la zone cible et son environnement. 

2. Ensemble de traitement thermique selon la revendication 1, caracterise par le 
fait que 1'unite de controle (300) comprend en outre des moyens (340) pour estimer les 
pertes locales en energie thermique, a partir d'une estimation de la conduction 
thermique et de la distribution spatiale de la temperature dans la zone cible et son 
environnement. 

3. Ensemble de traitement thermique selon Tune des revendications precedentes, 
-earacterise-par4e-fai^ emettent des ultrasons 

focalises. 

4. Ensemble de traitement thermique selon Tune des revendications precedentes 
caracterise par le fait que les moyens (200) pour mesurer et enregistrer la distribution 
spatiale de la temperature comprennent un appareil d'imagerie par resonance 
magnetique. 

5. Ensemble de traitement thermique selon Tune des revendications precedentes, 
caracterise en ce qu'il comprend des moyens devaluation de la distribution spatiale, 
dans la zone cible et son environnement, de l'energie fournie a la zone cible. 

6. Procede de regulation d'un ensemble de traitement thermique d'une zone cible 
d'un tissu biologique (410), comprenant I'etape consistant a deposer localement de 
l'energie dans la zone cible, 
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caracterise par le fait qu'il comprend en outre les etapes consistant a 

- evaluer les gradients de temperature dans la zone cible et son environnement ; et 

- en deduire Tenergie a deposer dans la zone cible pour atteindre la temperature 
souhaitee. 

7. Procede de regulation d'un ensemble de traitement thermique, selon la 
revendication 6, caracterise par le fait qu'il comprend en outre l'etape consistant a 
estimer les pertes locales en energie, dans la zone cible et son environnement. 

8. Procede de regulation d'un ensemble de traitement thermique, selon la 
revendication 7, caracterise par le fait qu'il comprend en outre une etape consistant a 
evaluer la distribution spatiale, dans la zone cible et son environnement, de Penergie 
fournie a la zone cible. 
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REVENDICATIONS 

1. Ensemble pour le traitement thermique d'une zone cible d'un tissu biologique 
(4 1 0), comprenant 

5 - des moyens generateurs d'energie (100) pour fournir localement de l'energie dans la 
zone cible ; 

- des moyens (200)pour mesurer et enregistrer la temperature dans la zone cible; 

- une unite de controle (300) comprenant des moyens (330) pour determiner, a partir de 
la temperature mesuree dans la zone cible, la quantite d'energie devant etre fournie a la 

10 zone cible, et des moyens pour commander (350) aux moyens generateurs de l'energie 
(100), de delivrer cette valeur de la puissance ; 

caracterise par le fait que l'unite de controle (300) comprend en outre des moyens de 
traitement numerique (320), point par point, de la distribution spatiale de la temperature 
dans la zone cible et son environnement. 
15 2. Ensemble de traitement thermique selon la revendication 1, caracterise par le 

fait que l'unite de controle (300) comprend en outre des moyens (340) pour estimer les 
pertes locales en energie thermique, a partir d'une estimation de la conduction 
thermique et de la distribution spatiale de la temperature dans la zone cible et son 
environnement. 

20 3. Ensemble de traitement thermique selon l'une des revendications precedentes, 

caracterise par le fait que les moyens generateurs d'energie (100) emettent des ultrasons 
focalises. 

4. Ensemble de traitement thermique selon Tune des revendications precedentes 
caracterise par le fait que les moyens (200) pour mesurer et enregistrer la distribution 

25 spatiale de la temperature comprennent un appareil d'imagerie par resonance 
magnetique. 

5. Ensemble de traitement thermique selon Tune des revendications precedentes, 
caracterise en ce qu'il comprend des moyens devaluation de la distribution spatiale, 
dans la zone cible et son environnement, de l'energie fournie a la zone cible. 
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